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摘 ”要 : 无 定 河 流域 是 黄土 高 原生 态 恢 复工 程 实施 的 重点 区 ,探究 其 土地 覆 被 的 空间 分 异 机 制 及 相关 水 碳 变量 的 
变化 特征 ,对 支撑 区 域 水 土 资源 保护 与 规划 以 及 服务 区 域 生态 文明 建设 等 工作 具有 重要 作用 。 本 文 使 用 线性 倾向 
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法 .Mann-Kendall 趋势 检验 Pettitt 突变 检验 HE FAS tie TT TS aT TC eT Yat Sat a BD 8 ECR IE, Be E E 


SK ZR , If WAH 7K he RB BS FE RT J EP (GPP) SEAS BA (ET) AKIA BK (WUE ) 等 关 
键 水 碳 变 量 的 变化 特征 。 研 究 表明 :(1) 1990 一 2019 年 期 间 ,流域 整体 草地 、 林 地、 建设 用 地 显著 增加 ,耕地 .荒地 显 
著 减 少 ,其 中 林 草 面积 增加 区 域 主要 集中 在 流域 下 游 及 无 定 河沿 岸 地 区 ; (2) 人 口 密度 .降水 .气温 等 对 流域 土地 覆 
被 空间 格局 具有 重要 影响 ,整体 而 言 社会 经 济 因素 的 影响 大 于 自然 因素 ,但 以 降水 气温 为 代表 的 自然 因素 的 影响 
在 增强 ;(3) 流域 水 碳 变 量 的 变化 与 土地 覆 被 变化 具有 和 较 好 的 对 应 关系 ,空间 上 ,以 耕地 林地、 草地 为 主要 覆 被 的 
流域 东南 部 的 GPP、.ET、WUE 相对 偏 高 ,以 草地 .荒地 为 主要 履 被 的 流域 西北 部 的 CPP ET WUE 相对 偏 低 , 时 间 上 ， 
2001 一 2019 年 间 , 流 域 整体 GPP、.ET、WUE 均 旦 增加 趋势 ,其 中 GPP、WUE 在 流域 绝 大 部 分 区 域 均 显著 增加 ,而 ET 


ul 


主要 在 流域 中 下 游 地 区 显著 增加 。 在 退耕 还 林 还 草 生态 恢复 工程 的 实施 与 气候 变化 背景 下 ,无 定 河流 域 林 草 得 到 


恢复 ,生态 环境 转 好 。 


关键 词 : 黄土 高 原 ; BOL; 实际 蒸 散 发 ;总 初级 生产 力 ; 水 分 利用 效率 ; 地 理 探测 天 


土地 歼 被 作为 连接 土壤 ` 大 气 水文 等 生态 要 
素 的 重要 纽 市 ,深刻 影响 陆 气 能 量 交换 .陆地 水 碳 
循环 和 区 域 生 态 环境 演变 。 作 为 全 球 环境 变化 的 
重要 组 成 部 分 ,土地 覆 被 变化 不 仅 能 反映 陆地 表 


早 半 干 早 区 ,是 我 国生 态 系统 最 脆弱 的 地 区 之 一 ， 
对 气候 变化 和 人 类 活动 十 分 敏感 ""。 其 中 无 定 河 
流域 作为 典型 的 生态 脆弱 地 区 "", 是 黄土 高 原 退 耕 
还 林 还 草 生态 恢复 工程 实施 的 重点 区 域 ,也 是 国家 


层 景观 的 时 空 动 态 过 程 ,还 可 以 在 不 同 程度 上 改变 
全 球 和 区 域 的 水 循环 能量 循环 和 生物 物质 循环 过 
程 ””。 因 此 ,国内 外 学 者 在 全 球 、 区 域 及 流域 尺度 
上 对 土地 履 被 动态 格局 .驱动 机 制 .生态 环境 效应 
等 方面 做 了 大 量 研 究 '" ,对 认识 气候 变化 方向 、 
人 -地 复杂 系统 以 及 地 -水 - 碳 相 互 作用 机 理 等 具有 
重要 参考 价值 。 

自 20 世 纪 70 年 代 以 来 ,中 国 在 黄土 高 原 地 区 
开展 了 大 规模 生态 保护 和 生态 恢复 工程 ,人 为 对 当 
地 土地 履 被 格局 进行 了 干预 六” 。 黄 土 高 原 位 于 于 
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贯彻 实施 黄河 流域 生态 保护 和 高 质量 发 展 战略 的 
重要 地 区 ,流域 内 人 -地 相互 作用 剧烈 , 覆 被 变化 明 
显 ,尤其 植被 恢复 显著 '”。 针 对 无 定 河流 域 以 及 整 
个 黄土 高 原 地 区 ,以 往 与 土地 覆 被 相关 的 研究 侧重 
分 析 土 地 覆 被 (例如 NDVI) 本 身 的 格局 变化 "“" ,或 
强调 覆 被 变化 在 改善 环境 水土 保持 等 方面 发 挥 的 
作用 和 导致 河川 径流 下 降 等 问题 “” ,而 考虑 地 - 
水 - 碳 耦 合 机 制 以 评估 生态 系统 关键 变量 的 研究 相 
对 薄弱 ;此 外 ,对 土地 和 覆 被 空间 分 异 机 制 的 定性 分 
析 多 ,定量 化 研究 及 机 理 研 究 相 对 不 足 。 因 此 , 综 
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合 评 价 土地 覆 被 的 时 空格 局 .空间 分 异 机 制 相关 
水 碳 变量 变化 特征 具有 重要 意义 。 

以 往 开展 土地 覆 被 时 空 变化 特征 分 析 时 ,针对 
长 时 间 连 续 数 据 序列 ,通常 采用 简单 线性 回归 Sen’ s 
斜率 估计 、Mann-Kendall 检 验 等 方法 描述 变量 的 时 
间 变 化 3; 针对 不 连续 数据 ,多 采用 转移 和 矩阵、 图 
谱 分 析 ,重心 迁移 模型 等 方法 分 析 土 地 履 被 的 转型 
结构 .时 空 演变 等 。 在 开展 土地 覆 被 变化 驱动 力 
分 析 时 ,通常 将 影响 要 素 划 分 为 自然 要 素 和 社会 经 
济 要 素 两 大 类 ,其 中 自然 要 素 影响 主要 表现 为 累积 
性 效应 ,而 社会 经 济 要 素 的 影响 效果 快速 而 显著 ”， 
以 往 研究 主要 通过 定性 描述 或 使 用 加 权 回 归 ,logis- 
tic 回归 模型 等 方法 对 驱动 力 进 行 探究 ”” ,然而 相 
关 人 研究 往往 主要 针对 土地 覆 被 的 时 间 变 化 ,对 其 空 
间 分 异 机 制 的 研究 相对 较 少 。 近 年 来 ,地 理 探测 需 
2 为 单 期 下 垫 面 数据 的 空间 驱动 机 制 分 析 提 供 了 
有 效 工具 ,在 量化 多 环境 因子 对 土地 覆 被 空间 格局 
的 影响 方面 获得 广泛 应 用 ”i。 在 水 碳 变量 分 析 方 
H , ZS AC (Evapotranspiration , ET) 连接 水 循环 和 能 
量 循环 等 关键 过 程 ,是 绿 水 资源 耗 散 状况 的 重要 指 
ERO, MW RÆ J (Gross Primary Productivity, 
GPP ) 72 26 (iE AG ESRAR COH EE ER, 
水 分 利用 效率 (Water Use Efficiency, WUE) ES 
生态 系统 碳水 关系 的 综合 指标 ”” ,这 3 个 水 碳 变 
量 作为 陆地 水 碳 循 环 的 关键 要 素 和 生态 系统 关键 
变量 ,对 土地 覆 被 变化 具有 敏感 响应 ,因此 常 被 用 
做 土地 覆 被 水 碳 效 应 的 分 析 指 标 *“” 。 然 而 ,由 于 
A H URE R ET 与 CPP 数据 会 导致 WUE 的 错误 佑 
算 ,加 上 在 数据 分 析 时 考虑 耦合 水 碳 过 程 的 遥感 产 
品 有 限 ,当前 研究 主要 侧重 对 单个 变量 的 分 析 , 
合 水 碳 变量 的 研究 较 少 。 

基于 此 ,本文 以 无 定 河流 域 为 研究 对 象 ,首先 
使 用 连续 的 土地 覆 被 数据 ,利用 Mann-Kendall 趋势 

仿 验 空间 分 析 等 方法 ,分 析 流 域 土地 覆 被 的 时 空 
变化 ,在 此 基础 上 ,结合 气象 .DEM、 人 口 .GDP 数 
据 ,运用 地 理 探测 带 探 究 多 个 自然 因子 和 社会 经 济 
因子 对 土地 覆 被 空间 分 异 的 影响 ,最 后 ,依靠 水 碳 
耦合 的 PML_V2 陆地 蒸 散发 与 总 初级 生产 力 数据 
集 ,在 同 源 数据 基础 上 分 析 流 域 蒸 散 发 .总 初级 生 
产 力 .水 分 利用 效率 的 时 空 变化 ,以 期 为 黄土 高 原 
地 区 的 水 土 资源 利用 、 生 态 文明 建设 等 提供 科学 
参考 。 


1 数据 与 方法 


1.1 研究 区 概况 

无 定 河 流域 位 于 37?09'~39”03'N 、107"45' ~ 
110°34’E Zh] , 跨 内 蒙古 .陕西 两 省 ,面积 约 为 3.6x 
10° km’, it Jat Ey FETE 834 ~ 2001 mm 之 间 ,西北 高 , 东 
南 低 , 其 地 理 位 置 及 高 程 如 图 1 所 示 , 西 北部 位 于 毛 
马 素 沙 地 南 缘 ,东南 部 位 于 黄土 高 原 丘 陵 沟 蜜 区 ， 
生态 脆弱 。 气 候 属 温带 大 陆 性 季风 气候 ,流域 多 年 
平均 气温 约 为 8.9 'C ,多 年 平均 降水 量 约 为 369.9 
mm, 年 内 分 配 不 均 , 主 要 集中 在 7 一 9 月 ,降水 空间 
分 布 为 东南 多 ,西北 少 ' 直 。 无 定 河 是 黄河 中 游 河口 
至 龙门 区 段 最 大 的 支流 ,发 源 于 陕西 省 榆林 市 定 边 
BAF WIGS MACY HA ONG Ae 
县 ,由 西北 向 东南 注入 黄河 ,干流 全 长 491 km, 主要 
支流 有 北岸 的 纳 林 河 ,海流 免 河 、 榆 溪 河 和 南岸 的 
Fm KET ETS 


108°E 109°E 


110°E 
图 例 


图 1 无 定 河流 域 地 理 位 置 及 高 程 


Fig. 1 Location and elevation ofthe Wuding River Basin 


1.2 数据 来 源 

土地 履 被 数据 (1990 一 2019 年 , 共 30 期 ) 来 自 中 
国土 地 覆盖 产品 ,空间 分 辩 率 为 30 m ,数据 基于 
5463 个 目 视 解 译 样 本 ,总 体 精度 达 80%。 本 文 主要 
关注 草地 、 耕 地、 林地 建设 用 地 .荒地 等 5 大 类 型 ， 
故 对 源 数据 进行 重 分 类 ,其 中 耕地 为 源 数据 中 的 农 
田 ,林地 为 源 数据 中 的 森林 和 灌木 ,建设 用 地 为 源 
数据 中 的 不 透水 区 ,草地 与 荒地 和 源 数 据 类 型 
相同 。 

年 平均 气温 年 降水 量 、 人 口 密 度 GDP ,高 程 匀 
数据 来 自 中 国 科学 院 资 源 环境 科学 与 数据 中 心 
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(https://www.resde.cn/) ,空间 分 辩 率 均 为 1 kmo Æ 
平均 气温 年 降水 量 、 人 口 密度 、GDP 等 数据 的 年 份 
为 1995 年 .2000 年 .2005 年 .2010 年 .2015 年 .2019 
年 ( 共 6 期 )”“” ,坡度 . 坡 向 数据 基于 高 程 数 据 , 通 
过 ArcGIS 软 件 计 算得 来 。 

蒸 散 发 (Evapotranspiration ,ET) 、 总 初级 生产 力 
(Gross Primary Productivity, GPP ) 数 据 来 自 “ 中 国 区 
Ja PML-V2 陆地 蒸 散 发 与 总 初级 生产 力 数据 集 
(2001 一 2019 年 , 共 19 期 )*”(http://data.tpdc.ac.cn/ 
zh-hans/), 数 据 空间 分 辨识 为 500 m, 基 于 包含 中 国 9 
种 植物 功能 类 型 的 26 个 涡 动 协 方差 通 量 塔 的 观测 
值 校准 显示 , ET 数据 的 Nash-Sutcliffe 效率 (NSE ) 为 
0.75, 均 方 根 误差 (RMSE) 为 0.69 mm: d7’; GPP žiji 
的 NSE 为 0.82,RMSE 为 1.71gC.m”…d'。 根据 定 
X ,水 分 利用 效率 (Water Use Efficiency, WUE) 为 
GPP 与 ET 之 比 。 

本 文 从 以 上 数据 集中 提取 了 无 定 河流 域 范 围 
内 的 数据 ,用 于 分 析 研 究 。 

1.3 研究 方法 

1.3.1 趋势 .突变 检 验方 法 ”本文 利用 线性 倾向 估 
计 、Mann-Kendall 趋势 检验 .Pettitt 突变 检验 来 分 析 
土地 履 被 及 水 碳 变量 序列 的 趋势 性 与 突变 点 。 其 
中 线性 倾向 估计 通过 求 变量 关于 时 间 的 线性 回归 
来 表征 变量 的 变化 趋势 ,其 中 回归 系数 可 以 反映 趋 
势 方 向 和 倾向 程度 ”“。Mann-Kendall 趋 势 检 验 不 需 
要 样本 数据 遵从 特定 的 分 布 ,能 够 有 效 的 别 除 异常 
值 ,适用 于 非 正 态 分布 数 据 的 趋势 检验 ””。 检 验 
统计 量 Z 表 征 数据 序列 的 趋势 性 ,121 值 越 大 ,趋势 越 
明显 ;Z > 0, 说 明 有 增 大 趋势 ,Z<0, 说 明 有 减 小 趋 
势 。 当 给 定 显 著 性 水 平 w 后 , 若 IZ 大 于 临界 值 , 则 序 
列 的 趋势 显著 ,本文 统一 选用 a = 0.05 的 显著 性 水 
平 ,所 对 应 区 间 为 +1.96。Pettitt 突变 检验 为 非 参 数 
检验 方法 ,最 初 由 Pettitt ”提出 ,用 于 检验 突变 点 。 
1.3.2 地 理 探测 器 ”地理 探测 器 ?是 探测 空间 分 异 
性 ,以 及 揭示 其 背后 驱动 力 的 一 组 统计 学 方法 , 包 
括 因子 探测 .交互 作用 探测 .风险 区 探测 和 生态 探 
测 4 部 分 。 其 中 因子 探测 可 用 于 探测 某 因子 X 对 属 
性 7 的 空间 分 异 的 解释 力度 ,该 力度 用 g 值 (0 < q< 
1) 度 量 ,4 值 越 大 表示 自 变 量 X 对 属性 了 的 解释 力 越 
强 , 计 算 公 式 如 下 : 


L 
SNo 
h=1 


No’ 
TU A= 1…… , 江 为 变量 了 或 因子 和 的 分 层 , 即 分 


q=1- 
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类 或 分 区 ;WV 和 NN 分 别 为 层 h 和 全 区 的 单元 数 ; oj 
Filo? 分 别 是 层 h 和 全 区 的 Y 值 的 方差 。 

本 文选 取 年 平均 气温 年 降水 量 、 坡 度 、 坡 向 、 
高 程 等 自然 因子 ,人 口 密度 .GDP 等 社会 经 济 因子 ， 
运用 地 理 探测 器 探究 不 同 自 然 .社会 经 济 因 子 对 土 
地 覆 被 空间 分 异 的 影响 。 首 先 , 利 用 最 优 离散 化 方 
TES? , 将 年 降水 量 、 年 平均 气温 ,高程 .坡度 D 
人 口 密度 、GDP 等 数值 型 数据 进行 离散 化 ,转化 为 
类 型 变量 ,其 次 ,利用 最 邻近 内 插 法 将 土地 覆 被 数 
据 的 空间 分 辩 率 重 采样 至 1 km ,最 后 利用 地 理 探测 
絮 对 土地 履 被 空间 分 异 驱动 进行 计算 。 


2 结果 与 分 析 


21 the RN Sey 

EE Ss AU AE Ss FT ETT it Jae EE AY) Sa BK 
类 型 ,1990 一 2019 年 期 间 ,流域 整体 草地 、 林 地 、 建 
设 用 地 显著 增加 ,耕地 荒地 显著 减少 , 林 草 面积 增 
加 区 域 主要 集中 在 流域 下 游 及 无 定 河沿 岸 地 区 , 建 
设 用 地 在 榆林 市 周边 明显 增加 ,耕地 主要 减少 在 流 
域 下 游 地 区 ,而 荒地 主要 减少 在 流域 中 上 游 地 区 。 
整体 看 来 ,1990 一 2019 年 ,无 定 河流 域 退 耕 还 林 还 
草 生 态 恢复 工程 成 效 显著 ,生态 环境 得 到 改善 。 

具体 而 言 ,如 图 2、 图 3 所 示 , 草 地 在 全 流域 内 广 
泛 分 布 , 多 年 平均 面积 占 流域 总 面积 的 68.62% , 草 
地 面积 在 1990 一 2019 年 间 呈 显著 增加 趋势 (Z = 
6.49,P< 0.05), 并 在 2002 年 存在 突变 点 ,增加 区 域 
主要 分 布 在 流域 下 游 地 区 以 及 沿 无 定 河水 系 附近 ; 
耕地 是 流域 内 面积 第 二 的 土地 覆 被 类 型 ,主要 分 布 
在 流域 中 下 游 , 多 年 平均 面积 占 比 为 19.39% , 虽然 
流域 中 上 游 局 部 地 区 出 现 新 增 耕 地 ,但 下 游 的 耕地 
大 面积 减少 ,导致 流域 整体 耕地 面积 在 1990 一 2019 
年 面积 显著 减少 (Z = -3.39,P < 0.05) ;荒地 面积 在 
流域 内 占 到 11.21% ,主要 分 布 在 流域 东北 部 ,在 
1990—2019 年 期 间 ,荒地 呈 显 著 减 少 趋势 (2Z = 
-6.42,P < 0.05) ,减少 地 区 主要 位 于 无 定 河 上 游 支 
流 河岸 附近 ,与 新 增 草地 区 域 有 较 好 的 一 致 性 ; 建 
设 用 地 面积 占 流域 面积 的 0.33% ,近年 来 面积 呈 增 
加 趋势 (Z = 7.74,P < 0.05), 尤 其 是 在 2004 年 之 后 迅 
速 增加 ,新 增 地 区 主要 集中 在 榆林 市 ;林地 在 流域 
内 占 比 较 少 ,为 0.03% ,但 呈 显 著 增 加 趋势 (2 = 
3.96,P < 0.05) ,尤其 在 2006 年 后 迅速 增加 ,新 增 林 
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图 2 1990 一 2019 年 无 定 河流 域 各 土地 覆 被 类 型 面积 变化 
Fig. 2 Area changes of land cover in the Wuding River Basin, 1990-2019 
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地 主要 集中 在 流域 的 最 南部 。 


2.2 土地 覆 被 空间 分 异 驱 动力 探测 
综 上 所 述 , 无 定 河 流域 各 土地 替 被 类 型 分 布 存 
在 明显 的 空间 分 异性 ,并且 不 同 区 域 土地 履 被 类 型 


108°E 109°E 110°E 


图 3 无 定 河流 域 1990 一 2019 年 土地 覆 被 类 型 变化 情况 
Fig.3 Changes of land cover in the Wuding River Basin, 1990 -2019 


还 存在 不 同方 向 的 转变 。 因 此 ,进一步 从 自然 、 社 


如 图 4 所 示 。 


会 经 济 等 多 角度 出 发 ,利用 地 理 探测 器 对 无 定 河流 
域 6 期 土地 莉 被 的 空间 分 布 驱动 进行 了 探测 , 结 
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图 4 土地 堆 被 空间 分 异 因 子 探测 结 
Fig. 4 Factor detection results of land cover 


spatial differentiation 


从 多 年 平均 解释 力 看 ,在 各 自然 因子 中 ,降水 
是 影响 无 定 河 流域 土地 覆 被 空间 分 布 机 制 的 主 控 
因素 ,在 社会 经 济 因 子 中 ,人 口 密度 比 GDP 的 影响 
力 大 。 有 具体 而 言 ,降水 .气温 .高程 .坡度 . 坡 向 对 无 
定 河流 域 土 地 覆 被 空间 分 异 的 多 年 平均 解释 力 g 值 
分 别 为 0.09 .0.08 .0.07 .0.06.0.02, 人 口 和 GDP 的 多 
年 平均 解释 力 g 值 分 别 为 0.12、0.07, 以 上 因子 对 人 研 
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究 区 土地 覆 被 空间 分 布 的 影响 均 通 过 显著 性 检验 。 

此 外 ,本 文 将 自然 .社会 经 济 方面 解释 力 最 大 
的 前 两 项 进行 加 和 ,得 到 分 别 代表 自然 和 社会 经 济 
因子 相对 解释 力 大 小 结果 的 时 间 变 化 曲线 。 由 图 4 
可 知 , 除 2000 年 .2005 年 外 ,社会 经 济 因 素 对 土地 者 
被 的 空间 分 异 解释 力 大 于 自然 因素 ,但 在 1995 一 
2019 年 间 ,社会 经 济 因素 的 解释 力 存在 微弱 下 降 趋 
势 ,而 自然 因素 的 解释 力 存在 微弱 上 升 趋势 , 相 比 
1995 年 ,二 者 间 的 差 中 在 减 小 。 
2.3 相关 水 碳 变量 的 时 空 变化 

在 自然 与 社会 经 济 等 多 方面 因素 的 影响 下 ,无 
定 河流 域 的 土地 覆 被 分 布 情况 在 2000 年 后 发 生 了 
明显 的 改变 ,这 在 一 定 程度 上 势必 会 影响 到 无 定 河 
流域 水 碳 变量 的 时 空 变化 。 因 此 ,本 文选 取 GPP、 
ET WUE 这 3 个 重要 水 碳 变量 展开 分 析 ,结果 如 图 5 
所 示 。 

从 空间 分 布 上 看 ,如 图 Sa BI Sb. BI Sc 所 示 ,无 
定 河流 域 多 年 平均 CPP ET, WUE 的 空间 分 布 与 土 
地 覆 被 空间 分 异 有 较 好 的 对 应 性 ,大 体 存在 由 东南 
向 西北 递 变 的 规律 。 东 南部 以 耕地 、 林 地 、 草 地 为 


(a) 多 年 平均 GPP (b) 多 年 平均 ET (c) 多 年 平均 WUE 
108°E 109°E 110°E 108°E 109°E 110°E 108°E 109°E 110°E 
z z Z T T T Zz Z, ii T T Z 
A a S BS a 
Z Zz g E g A 
加 B A ag S 
Z Ze > 0 30km | £ Z Z 
108°E 109°E 110°E 108°E 109°E 110°E 108°E 109°E 110°E 
(d) GPP 趋 势 变化 (e) ET 趋势 变化 (f) WUE 趋 势 变化 
108°E 109°E 110°E 108°E 109°E 110°E 108°E 109°E 110°E 
z, ZZ i ZZ Zz 
a -| a & 上 各 
foal nm mom M 
Z Z zZ Zz z 
A | % % œ oo f 20 
M m r wi i nn 
rs oS a ee 
z Zz Zz jans Z 
人 mt 
一 1.96 0 1.96 —1.96 0 1.96 一 1.96 0 1.96 
| 1 | 1 
108°E 109°E 110°E 108°E 109°E 110°E 108°E 109°E 110°E 


图 5 2001 一 2019 年 无 定 河流 域 GPP、ET、WUE 的 变化 
Fig. 5 Changes of GPP, ET and WUE in the Wuding River Basin, 2001-2019 
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图 6 2001—2019 年 无 定 河流 域 整体 及 不 同 土地 履 被 类 型 的 GPP .ET、WUE 变化 
Fig.6 Changes in GPP, ET and WUE of the whole and different land cover in the Wuding River Basin, 2001-2019 


EA ok HH KI GPP ET WUE 相对 偏 高 ,西北 部 
以 草地 .车 地 为 主要 履 被 的 地 区 其 CPP ET、WUE 相 
对 偏 低 。 此 外 ,水 系 附近 的 CPP、ET、WUE 相对 
偏 高 。 

从 时 间 变 化 上 看 ,如 图 6、 图 5e、 图 5f、 图 5g 所 
ZR ,2001 一 2019 年 间 ,流域 整体 CPP ET、WUE 均 呈 
增加 趋势 ,CPP .WUE 在 流域 绝 大 部 分 区 域 均 显著 
增加 ,而 ET 主要 在 流域 中 下 游 地 区 显著 增加 。 进 一 
步 分 析 不 同 土地 履 被 类 型 的 水 磋 变 量 及 其 变化 可 
以 发 现 ,林地 、 耕 地、 草地、 荒地 等 土地 履 被 类 型 的 
GPP, ET, WUE 也 且 增 长 趋势 ,林地 的 GPP, ET, 
WUE 在 各 土地 覆 被 类 型 中 占据 第 一 ,耕地 、 草 地 次 
之 ,荒地 最 小 。 此 外 ,不 同 土地 覆 被 类 型 的 CPP、 
WUE 差 别 较 大 ,并 且 不 同 土地 履 被 类 型 间 差 距 还 在 
扩大 ,而 不 同 土地 覆 被 类 型 间 的 ET 差别 较 小 。 


3 讨论 


无 定 河流 域 整体 CPP .ET、WUE 均 呈 增加 趋势 ， 
主要 受到 人 为 改变 土地 覆 被 与 气候 变化 的 影响 。 
植被 可 通过 光合 作用 吸收 降低 空气 中 C0, 浓 度 ,与 
其 他 土地 利用 类 型 相 比 ,植被 覆盖 较 大 的 林 草 地 的 


固 碳 能 力也 相对 较 强 。1999 年 ,陕西 作为 全 国 试 
点 ,首先 开展 了 退耕 还 林 还 草 生态 恢复 工程 ,故此 
大 大 促进 了 当地 GPP 的 增长 ,这 类 碳 固 存 增加 主要 
体现 林 草 地 恢复 明显 的 流域 南部 地 区 。 而 全 流域 
大 范围 CPP 增长 可 能 由 气候 变化 驱动 ,依据 主要 有 
两 点 ,一 是 空间 驱动 力 分 析 结 果 显 示 降 水 和 气温 对 
流域 土地 覆 被 格局 的 影响 在 逐渐 扩大 ,毛乌素 沙 地 
的 暧 湿 化 “在 一 定 程度 上 驱动 了 流域 西北 部 低 GPP 
的 荒地 向 高 CPP 的 草地 的 转变 ;二 是 不 同 土地 履 被 
类 型 的 GPP 均 在 增加 ,在 气候 变 暧 推动 下 所 有 土地 
窗 被 类 型 中 绿色 植物 的 生存 环境 得 到 改善 , 碳 固 存 
能 力 均 得 到 提高 。 与 GPP 相 比 ,流域 内 ET 显著 增加 
的 区 域 较 少 , 主要 分 布 在 中 下 游 地 区 , 而且 不 同 土 
地 禾 被 类 型 间 的 ET 差别 小 ,这 与 流域 的 自然 气候 条 
件 密切 相关 ,无 定 河流 域 位 于 干旱 半 干 旱 区 ,当地 
的 实际 蔡 散 发 量 受 到 水 分 供应 的 限制 *, 此 外 ,水 
分 还 是 制约 植被 恢复 的 主要 因素 。 因 此 ,水 分 条 件 
相对 优 瀑 的 流域 中 下 游 地 区 植被 恢复 明显 ,同时 
MEHE T AKRE. WUE A T GPP E 
ET 的 变化 特征 ,在 GPP 与 ET 均 以 增加 为 主 的 情 
下 ,流域 绝 大 多 数 区 域 的 WUE 仍 呈 显 著 增加 趋势 ， 
说 明 该 区 域 CPP 的 显著 增加 主导 了 WUE 的 增加 ,在 
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植被 恢复 耗 散 水 分 增加 的 情况 下 ,流域 整体 的 水 分 
利用 效率 是 提高 的 。 

在 无 定 河流 域 所 在 的 水 资源 不 足 的 黄土 高 原 
区 以 及 广大 西北 干旱 半 干 星 区 开展 生态 保护 工程 
实践 时 ,需要 以 水 资源 植被 承载 力 为 依据 ,权衡 植 
被 恢复 的 碳 效 益 和 水 效应 。 从 生态 环境 角度 看 , 植 
树 造林 会 增加 区 域 固 碳 能 力 ,改善 生态 环境 ;而 从 
水 资源 耗 散 角度 看 ,森林 履 被 增加 会 导致 更 强 的 蒸 
散发 ,水 资源 耗 散 显著 增加 ,根据 水 量 守恒 定律 , 那 
么 人 类 赖 以 利用 的 径流 量 就 会 减少 。 因 此 ,耦合 水 
碳 机 制 来 对 区 域 植树 造林 的 可 行 性 和 退耕 退 林 还 
草 效果 的 可 持续 性 进行 科学 评 佑 具有 重要 意义 。 
在 追求 林 草 植被 恢复 的 固 碳 功能 的 同时 ,有 必要 充 
分 考量 当地 水 分 条 件 .水 分 情势 变化 和 未 来 可 持续 
性 等 ,真正 做 到 以 水 定 绿 、 以 水 定 碳 , 因 地 因 水 开展 
适宜 性 生态 建设 工作 , 切 不 可 为 了 短期 内 提升 碳 汇 
能 力 而 忽视 水 分 条 件 的 限制 。 

本 文 在 无 定 河流 域 开展 土地 覆 被 的 变化 特征 
分 析 、 空 间 分 异 因 子 探测 、 相 关 水 矶 变量 变化 分 析 
的 “三 段 式 "综合 研究 模式 可 为 黄土 高 原 区 、 黄 河流 
域 及 其 他 区 域 的 地 -水 - 碳 耦 合 人 研究 提供 科学 参考 ， 
对 文 撑 区 域 水 土 资 源 保护 与 规划 实施 水 土 保持 工 
程 . 辅 助 实现 “ 双 碳 "目标 等 工作 也 具有 重要 的 实践 
参考 意义 。 但 本 人 研究 存在 土地 徐 被 空间 分 异 机 制 
分 析 不 连续 、 相 关 水 碳 变 量 分 析 数 量 不 全 面 等 不 
足 , 今 后 将 进一步 获取 更 高 精度 .逐年 连续 的 社会 
经 济 数据 .拓展 更 多 生态 环境 变量 等 来 展开 更 深入 
的 研究 。 


本 文通 过 分 析 1990s 以 来 无 定 河流 域 土地 覆 被 
的 基本 特征 、 空 间 分 异 驱 动力 以 及 相关 水 碳 变 量 
(GPP .ET、WUE ) 的 变化 特征 ,得 出 以 下 结论 : 

(1) 近 30 a 来 流域 林 草 覆 被 得 到 良好 恢复 , 退 
耕 还 林 还 草 生 态 恢 复工 程 的 实施 以 及 气候 变化 均 
是 其 背后 重要 的 推力 。 作 为 流域 内 最 主要 的 土地 
履 被 类 型 ,草地 于 2002 年 左右 出 现 明 显 的 增长 ,在 
流域 中 下 游 及 无 定 河 沿岸 地 区 大 面积 增加 ,其 中 下 
游 地 区 草地 主要 由 耕地 转化 而 来 ,中 上 游 河 岸 附近 
的 草地 主要 由 荒地 恢复 而 来 ;而 林地 明显 增加 出 现 
在 2006 年 左右 ,增加 地 区 也 仅 局 限 在 流域 最 南部 。 
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(2) 流域 土地 覆 被 在 空间 上 主要 呈现 由 东南 到 
西北 逐渐 变化 的 分 布 特征 ,人 口 密度 和 降水 对 该 格 
局 形成 的 影响 力 最 大 ,GDP 气温 ,高程 .坡度 . 坡 向 
等 因素 也 具有 重要 作用 。 此 外 ,降水 .气温 等 变 率 
较 大 的 自然 因素 对 流域 土地 覆 被 空间 分 异 的 影响 
旺 增强 趋势 。 

(3) 流域 GCPP、ET、WUE 的 空间 分 布 与 土地 覆 
被 空间 分 异 有 较 好 的 对 应 性 ,并 且 在 林 草 恢复 为 主 
的 覆 被 变化 背景 下 ,流域 整体 的 CPP ET WUE 明显 
增加 。 在 流域 绝 大 部 分 区 域 碳 固 存 增加 、 水 分 利用 
效率 提高 的 同时 ,伴随 着 中 下 游 地 区 水 分 耗 散 加 剧 
的 问题 ,但 总 体 而 言 ,流域 生态 环境 偏向 转 好 。 
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Spatial differentiation mechanism of land cover and related changes in 


water-carbon variables in Wuding River Basin 


LYU Jinxin'’*, LIANG Kang', LIU Changming', ZHANG Yihui'’’, LIU Lu 
(1. Key Laboratory of Water Cycle and Related Land Surface Processes, Institute of Geographic Sciences and Natural 
Resources Research, Chinese Academy of Sciences, Beijing 100101, China; 2. University of Chinese Academy of 


Sciences, Beijing 100049, China) 


Abstract: The Wuding River Basin is the key area for the implementation of the grain for green program on the 
Loess Plateau. Exploring the spatial differentiation mechanism of land cover and the characteristics of variation 
of related water-carbon variables is essential for supporting water and soil resources conservation and planning, 
along with serving the construction of regional ecological civilization. Here we used the linear tendency method, 
Mann-Kendall trend test, Pettitt test, and geodetector to analyze the spatiotemporal characteristics and the factors 
driving the spatial differentiation of land cover in Wudinghe River Basin. In addition, from the perspective of 
land- water- carbon coupling, we analyzed the characteristics of variation of total primary productivity (GPP), 
actual evapotranspiration (ET), and water use efficiency (WUE). Three main results were as follows: (1) The total 
grassland, forest, and construction land increased significantly, while the cropland and barren land decreased 
significantly during 1990-2019. The area of forest and grassland increased mainly in the lower reaches of the 
basin and along Wuding River. (ii) Population density, precipitation, and temperature have a significant impact on 
the spatial pattern of land cover. On the whole, the influence of socioeconomic factors is greater than that of 
natural factors, but the influence of natural factors represented by precipitation and temperature is increasing. (iii) 
There is a good correlation between the variation of water-carbon variables and land cover change. Spatially, 
GPP, ET, and WUE are relatively high in the southeastern part of the basin where cropland, forest, and grassland 
are the main cover, while they are lower in the northwestern part of the basin where grassland and barren land are 
the main cover. In terms of time, GPP, ET, and WUE all showed increasing trends during 2001-2019. GPP and 
WUE increased significantly in most parts of the basin, while ET increased significantly mainly in its middle and 
lower reaches. In conclusion, under the implementation of the grain for green program and climate change, the 
forest and grass in Wudinghe River Basin have been restored and the ecology has been improved. 

Keywords: the Loess Plateau; land cover change; actual evapotranspiration; gross primary productivity; water 
use efficiency; geodetector 


